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ABSTRACT 

We investigate the effect of photodesorption on the snow line position at the surface of a proto- 
planetary disk around a Herbig Ae/Be star, motivated by the detection of water ice particles at the 
surface of the disk around HD142527 by Honda et al. For this aim, we obtain the density and tem- 
perature structure in the disk with a 1+1D radiative transfer and determine the distribution of water 
ice particles in the disk by the balance between condensation, sublimation, and photodesorption. We 
find that photodesorption induced by the far-ultraviolet radiation from the central star depresses the 
ice-condensation front toward the mid-plane and pushes the surface snow line outward significantly 
when the stellar effective temperature exceeds a certain critical value. This critical effective temper- 
ature depends on the stellar luminosity and mass, the water abundance in the disk, and the yield 
of photodesorption. We present an approximate analytic formula for the critical temperature. We 
separate Herbig Ae/Be stars into two groups on the HR diagram according to the critical temperature; 
one is the disks where photodesorption is effective and from which we may not find ice particles at the 
surface, and the other is the disks where photodesorption is not effective. We estimate the snow line 
position at the surface of the disk around HD142527 to be 100-300 AU which is consistent with the 
water ice detection at > 140 AU in the disk. All results depend on the dust grain size by a complex 
way and this point requires more work in future. 

Subject headings: methods: numerical — protoplanetary disks — radiative transfer — stars: individual 
(HD142527) 



1. INTRODUCTION 

In a protoplanetary disk, there is the so-called snow 
line which is defined a s the inner bo undary of the ice- 
condensing region (e.g..[Ha"vashi 1981). Beyond the snow 
line, the surface mass density of solid materials is en- 
hanced due to condensation of water ice, and hence, the 
formation efficiency of large gaseous and icy planets is 
expected to be enhanced. Thus, the snow line plays a 
very important role in the planet formation. 

In spite of its importance, no one has detected the 
snow line position in a circumstellar disk so far. This is 
because the snow line is usually located close to the cen- 
tral star (R < 10 AU), where the spacial resolving power 
of current facilities is limited, and the disk is optically 
too thick to see through. On the other hand, the pres- 
ence of water ic e in circumstellar disks have already been 
detected (e.g.. lMalfa.it et al.l Il999t ITerada et all 120071: 
iHonda et all 120091 ). iTerada et all (|2007D observed scat- 
tered radiation from an edge- on disk and detected an 
absorption feature of water ice at 3 /im in it. However, 
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their observations give no information ab out the posi- 
tion o f the snow line along the line of sight. IHonda et all 
(2009) observed scattered light from the surface of an 
face-on disk and detected the 3 /im absorption feature 
by icy dust particles at R > 140 AU. They could es- 
timate the dust size as ~ 1 Mm by us i ng a c olor dia- 
gram metho d prop osed by llnoue et all (|2008f) . While 
IHonda et all (|2009b could not detect the snow line posi- 
tion, their method holds a potential to achieve it due to 
a high spatial resolution of their detection technique. 

Although an important factor for the planet formation 
is the snow line position at the mid-plane, detecting the 
snow line at the disk surface will be a first step because 
the disk surface is easily observed. The surface snow line 
is connected to the mid-plane snow line t hrough a 2D 
structure of the ice-condens ing region (e.g.. lDavyl2005l : 
Oka. Nakamo to fc Iddl27)Tl . If one knows the mecha- 
nism which determines the 2D structure of the snow line, 
one may convert the surface snow line into the mid-plane 
snow line. 

The main difference between the interior and the sur- 
face of protoplanetary disks is that the latter is exposed 
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to ultra-violet (UV) radiation from the central star and 
the interstellar radiation field. This UV radiation may 
affect the distribution of water ice at the disk surface. 
Indeed, it has been proposed that a non-thermal desorp- 
tion process is necessary for explaining the observed high 
abundance of gas-phase water molecules in cold outer re- 
gions of circumstellar disks and molecular clouds (e.g., 
iDominik et"aIll2Q05l:lHollenbach et~aTll2009H . A probable 
mechanism responsible for such a non-thermal desorption 
of water molecules is the photodesorption (or photosput- 
tering) process caused by far-UV (FUV) radiation. The 
photodesorption process has bee n recently studied the o- 
retically and experimentally fe.g..lWestlev et al.lll995aHb ; 
lAndersson et al.ll2006t lAnder sson fe van Di shocck 2008; 
lOberg et al.1l2009t lHama et al.l[2010h . 

Oberg et al. (2009) experimentally studied the pho- 
todesorption process and derived the photodesorption 
yield as a function of ice temperature and thickness of 
the ice layer. According to their results, the photodes- 
orption yield with ice temperature T and thickness of 
the ice layer larger than 8 mono layers is V p d (T, x > 
8) = 10" 3 (1.3 + 0.032T) UV photon" 1 , where x is the 
ice thickness scaled by a mono layer. They then car- 
ried out a numerical simulation about the distribution of 
gas phase H2O in a protoplanetary disk around a Herbig 
Ae/Be star, taking into account adsorption and thermal 
desorption processes on dust grains, and the photodes- 
orption process, and found that photodesorption reduced 
the H 2 ice abundance in the disk surface by orders of 
magnitudes. 

iGrigorieva et al.l (|2007t ) also showed theoretically that 
the presence of icy dust particles in a debris disk around 
a Herbig Ae/Be star is affected by the photodesorption 
process. The environment of the s urface of protoplane - 
tary disks such as one observed by iHonda et all (]2009) 
may be similar to that of debris disks, and hence the po- 
sition of the surface snow line would be affected by the 
photodesorption process. 

Recently, observations of water vapor in protoplane- 
tary disks show a rapid progress with Spitzer Space Tele- 
scope and Herschel Space Obs ervatory as well as large 
ground- based telescopes fe.g.. iPontoppidan et al.l 120101: 
Sturm et al.ll2010t Ivan Kempen et al.ll2010l : iBergin et al.l 
2010t iFedele et al.l 1201 11 1. There is lack of hot and 
warm vapor traced by near- and mid-infrared obser- 
vations in disks around Herbig Ae/Be stars, whereas 
such yapor was often found in disks around T T auri 
stars (|Pontoppidan et al"1l2010t IFedele et al.l 120111) . A 
difference of photochemistry is suggested as a cause of 
this difference. In Herbig Ae/Be disks, H2O molecules 
(photo) desorbed from icy grains may be photodisso- 
ciated by st r ong U V ra diation. On the o t her h and, 
iSturm et all (|2010fl and Ivan Kempen et al.l (|2010D de- 
tected cool water vapor in far-infrared observations in 
Herbig Ae/Be disks, suggesting incomplete photodisso- 
ciation in the outer part of the disks. More interestingly, 
IBergin et al.l (|2010f ) reported no detection of cold water 
vapor in submillimetre in a T Tauri disk and suggested 
that photodesorption was less efficient than theoretically 
expected. In any case, more works about photodesorp- 
tion and photodissociation in protoplanetary disks are 
required to depict a full picture of the photochemical 
evolution in the disks. In particular, less efficient pho- 



todesorption suggested by IBergin et al.l (|2010f ) should be 
tested also in Herbig Ae/Be stars. 

In this study, we numerically investigate the effect of 
the photodesorption process induced by the FUV radi- 
ation from the central star and the interstellar radia- 
tion field on the surface snow line position in circum- 
stellar disks. We focus on Herbig Ae/Be disks where 
photodesorption is effective as expected, while the pa- 
rameter space explored in this paper also covers T Tauri 
disks. The goal is to uncover what the determining pro- 
cess of the snow line position at a disk surface is. We 
omit photodissociation in this paper, but in fact, we find 
that this process is much slower than photodesorption 
(see §4.3.1). 

This paper is organized as follows. In §2, we describe 
our model for obtaining the distribution of icy dust par- 
ticles in a Herbig Ae/Be disk with photodesorption. In 
§3, we show effects of photodesorption on the position of 
the surface snow line and the dependence of the results 
on the stellar parameters. In §4, we consider an impli- 
cation of our results for observations, derive an analytic 
formula of the critical stellar effective temperature for 
an efficient photodesorption, and discuss effects of some 
processes omitted in our calculations. Finally, we sum- 
marize this paper in §5. 

2. MODEL 
2.1. Disk structure 

We consider a flared but geometrically-thin and 
optically-thick disk revolving around a central star in Ke- 
pler velocity. We assume that the disk is axi-symmetric 
about the rotational axis and plane-symmetric about the 
mid-plane. We adopt a coordinate system (R,Z), where 
R and Z are the radial distance from the central star and 
the vertical height from the disk mid-plane, respectively. 
The central star sits on the origin of the coordinate. 

We assume that the radial distribution of the surface 
density of th e disk is a power law distribution (e.g., 
IHavashi|[l98l as 

m) = siau \ a) 

where Siau is the surface density at R = 1AU and p is 
the power-law index describing the radial density distri- 
bution. We as sume Siau = 1700 g cm" 2 and p = 1.5 
(jHayashil 119811 ) throughout this paper. The disk is as- 
sumed to be in the hydrostatic equilibrium vertically, 
and heated by the irradiation by the central star. The 
temperature and density structures of the disk are nu- 
merically sol ved by the 1+ ID r adiati ve transfer method 
described bv iDullemond et ahl (|2002l ). We assume that 
tem peratures of dust and gas are always same. Un- 
like IDullemond et~a l. (2002), we consider isotropic scat- 
te ring of both stella r an d diffuse radiations accordin g 
to llnoue et all (|2009D and lOka. Nakamoto fc Idal (|2011l) . 
The details of the calculation are describ ed in Appendix 
A (see also lOka. Nakamoto fe Idd[2^TTh . Note that we 
assume that the grazing angle (j3) of the stellar radiation 
incident on the disk surface is uniformly 0.05 throughout 
the disk for simplicity. This is fairly artificial but it was 
very difficult to determine both the grazing angle and 
the ice abundance simultaneously in our 1+1D disk with 
scattering by an iteration scheme because of very slow 
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Figure 1. Absorption and scattering coefficients of the disk gas 
containing 1 /^m-sized silicate and water ice particles adopted in 
our calculation. These coefficients arc taken from Miyake & Naka- 
gawa (1993). Note that all the water molecules are contained in the 
icy dust particles. The thick and thin curves represent the coeffi- 
cients by icy and silicate dust particles, respectively. The solid and 
dashed curves represent the absorption and scattering coefficients, 
respectively. 

convergence. The mean molecular weight of the disk gas 
is assumed to be 2.3. 

Although we fix global parameters of the disk such as 
Siau, P, and /3, it would not affect the results presented 
in this paper very much when discussing the disk sur- 
face where the radiative transfer effect is negligible and 
the physical condition is determined almost locally. This 
enables us to derive an analytic formula describing the 
numerical results in §4.2, while the absolute scaling may 
be affected slightly. 

2.2. Dust opacity 

A very important factor for determining the structure 
and the radiation field in the disk is the optical property 
of the disk. This is mainly determined by the absorp- 
tion and scattering processes by dust particles. The gas 
opacity is negligible for the radiation considered in this 
paper, except for Lya photons. We omit gas opacity 
throughout this paper for simplicity, whereas an impor- 
tance of Lyq transfer was sugg ested (iFogel et al.1 120111 : 
iBethell fc~ Bcrgin 2011). It is still reasonable when dis- 
cussing the disk surface where the transfer effect is small. 

We consider only pure crystalline silicate and pure 
crystalline water ice particles floating in the disk gas 
as dust. These particles are assumed to be spherical 
with a fixed uniform size to avoid complexity. Here we 
assume 1 /xm-sized particle s as found in the HD142527 
disk by iHonda et al.l ([20095, but the effect of the grain 
size is discussed in §4.3.4. The mass fractions of sili- 
cate and water to the dis k gas are assumed to be 0. 0043 
and 0.0094, respectively, (jMivake fc Nakagawa! 1993h and 
uniform throughout the disk. This water abundance cor- 
responds to the number fraction of H2O molecules rela- 



tive to all gas particles in the disk as X 



h 2 o 

total 



1.2 x 10" 



which is the maximum fraction of water vapor Q The 
absorption and scattering coefficie nts of silicate and wa- 
ter ice particles are also taken from lMivake &: Nakagawal 
(1993). The coefficients for 1 /xm-sized particles for a 



The H2O abundance in Mivake & Nakagawa (1993) was not 
taken into account th e recent downward revision of O abundance 
l|Asplun d ct al. 2009T). 



unit disk mass are shown in Figure 1. 

2.3. Ice abundance 

We evaluate the abundance of water ice in the disk 
by the balance between changing rates of the ice dust 
radius due to condensation s con , sublimation s su b, and 
photodesorption s p( j; 



Ssuh(T) + Scon(r, A v " OI . 



s P d = 0, 



(2) 



where X^^ r is the number fraction of water vapor of 
the disk medium and T is the dust (and gas) temper- 
ature. The rates are evaluated by the formulas from 
iGrigorieva et al.l (|2007l ) (see Appendix B). The dust tem- 
perature T is obtained by the radiative transfer calcula- 
tion described in §2.1 and Appendix A. 

The number fraction of water vapor X^® T is deter- 
mined as equation ([2]) is fulfilled. Then, the number 

TT Q 

fraction of H2O in water ice to the disk gas A ic 2 is 
given by the residual of the subtraction of Xy^ T from 
the total number fraction of H2O A^^J. On the other 
hand, if s sub + s con + s pd < with A^ 2 p °. = A t " 2 °, wa- 
ter ice does not condense and X^° — 0. Using A t ^ 2 ^ 

and Xfc2 , the mass ratio of water ice to the total water 
is given by 



ice 

abundance x 



X 



h 2 



X 



H 2 ' 

total 



(3) 



Then, the absorption and scattering coefficients of a unit 
disk mass containing dust particles, K aDS and 

^sca are 

given by 



^abs — ^sil,abs ~t~ ^icc^icc,abs; (4) 
^sca — ^sil,sca ~t~ ^icc^icc,sca; (5) 

where K s ii,abs ; K sii,sca ; K ice,absj and Ki C e,sca are the ab- 
sorption and scattering coefficients of silicate and wa- 
ter ice dust particles bv lMivake fc Na kagawa (1993|), re- 
spectively (Kice.abs and Ki COjSca are the coefficients in the 
case where all the water molecules are condensed as icy 
dust particles). Since x- lcc itself affects the disk struc- 
ture through the radiative transfer, we solve x- lce and the 
disk structure iteratively to obtain consistent solutions 
of them. 

2.4. Stellar and interstellar FUV radiation 

We consider three sources of UV radiation: the stellar 
photosphere, an excess stellar radiation (probably caused 
by an accretion shock), and the interstellar radiation 
field. The stellar photospheric radiation is treated in the 
radiative transfer code described in §2.1 and Appendix 
A. We assume the spectrum to be a Planck function. The 
transfer of the other two radiation sources are treated as 
described in the subsequent paragraphs. As found in §2.5 
and §3.2, however, the two sources have a very weak ef- 
fect on the position of the surface snow line. The heating 
by the two sources is taken into account in equation (A8) 
in addition to that of the photospheric radiation. 

An excess FUV radiation is often observed in T Tauri 
stars, while the p hotospheric radiati on is dominant in 
Herbig stars (e.g.. iValenti et al1l2003[ ). Nevertheless we 
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consider the excess radiation for discussions of the de- 
pendence of the stellar effective temperature in §3.2. We 
assume the excess FUV luminosity to be a fraction of 
to the bolometric luminosit y of the photospher e 
(L*) as observed in T Tauri stars (jHerczeg et al.l [2004). 
This reproduces the observational trend that the ex- 
cess relative to the photospheric radiation is less impor- 
tant for earlier spectral type. Observations show that 
the excess FUV radi ation is dominated by Lya photons 
(jHerczeg et al.l 120041 ) . We assume that the excess radi- 
ation consists only of Lya photons for simplicity. How- 
ever, we omit any scattering of the Lya photons (see §2.2) 
although the resonant scattering of the photons is impor- 
tant in the sub-surface layer of the disk ([Fogel et"aLll2011fc 
iBethell fc^ Bcrgin 2011). The Lya assumption is only af- 
fects the calculation of the optical depth in the disk and 
the photon number flux. Consequently, the excess FUV 
radiation flux in the disk can be obtained as 

.excess , excess, 

(6) 

where i<FUV, excess, o{R) = 10 3 L*/(47ri? 2 ) is the incident 
excess FUV flux from the central star at the radius R for 
no dust case, tftjy(Z) is the absorption optical depth at 
the FUV (Lya) wavelength measured along the Z-axis 
from (R, +oo) to (R, Z), and (3 is the grazing angle for 
the stellar radiation. 

We also consider the FUV radiation from the interstel- 
lar field . We a dopt an average FUV radiation flux from 
iHabind jT968). The FUV flux is evaluated by assuming 
that it penetrates into the disk parallel to the Z-axis: 

F F tjvas{R,Z) = F FUV ,is,oexp(-T FU v(2')), (7) 

where -FFUV,is,o(-Rj Z) is the incident interstellar FUV ra- 
diation flux and is set to be 1.6 x 10 -3 erg cm -2 s _1 . The 
interstellar FUV radiation is also assumed to be con- 
sisted only of Lya photons as the stellar excess FUV. 
This is simply in order to treat the two additional FUV 
radiation sources homogeneously, whereas the interstel- 
lar radiation may be described by a Planck function or 
a power-law function made by a combination of Planck 
functions with various temperatures. 

2.5. Examples of disk vertical structure 

Figure [2] shows examples of the vertical structure of 
two annuli at R — 10 AU (left-hand tandem panels) and 
100 AU (right-hand tandem panels) around a central star 
with the effective temperature T* = 5000 K, the luminos- 
ity L* = 10 Lq, and the mass M* = 1 M . The top pan- 
els show the temperature structure along the disk vertical 
axis. The middle panels show the changing rates of the 
grain size by condensation (double-dot-dashed line), sub- 
limation (dotted line), and photodesorption (solid line). 
We also show the changing rates by photodesorption due 
to stellar (photospheric and excess) radiation (dashed 
line) and due to interstellar radiation (dot-dashed line). 
The bottom panels show the ice condensation fraction. 

In both annuli, condensation is rapid near the mid- 
plane and virtually all water is in ice. As the tem- 
perature increases (and the density decreases) along the 
vertical height, condensation becomes slower and slower. 
On the other hand, sublimation and photodesorption be- 
come suddenly effective beyond a certain height. Then, 
the ice condensation front is formed at a height where 



condensation balances with sublimation or photodesorp- 
tion. In the R = 10 AU case, sublimation is faster than 
photodesorption at the front, and then, the balance is es- 
tablished between condensation and sublimation. On the 
other hand, in the R = 100 AU case, sublimation is more 
than ten orders of magnitudes slower than photodesorp- 
tion because of a low temperature. Thus, condensation 
balances with photodesorption. 

Photodesorption by stellar FUV radiation always dom- 
inates that by interstellar FUV radiation at the disk sur- 
face because of a much stronger flux of the former ra- 
diation. However, the latter dominates the former at a 
sub-surface layer. This is caused by the different incident 
angles of the two radiations; the stellar radiation enters 
the disk with a small angle against the surface (i.e. graz- 
ing angle f3) and the penetration depth becomes shallow, 
whereas the interstellar radiation enters the disk verti- 
cally and reaches a deeper layer. When discussing the 
surface snow line, we conclude that the interstellar ra- 
diation hardly contributes to the determination of the 
surface snow line. Therefore, we omit the contribution 
by the interstellar FUV radiation in the remaining part 
of this paper. 

2.6. Definition of the surface snow line 

Let us define the position of the surface snow line, 
here. We consider the method observing the surface 
snow line by using the 3 /im H2O ice absorption fea- 
ture in the scattered light propos ed bvllnoue et al.l (120081 ) 
theoretically and attempted bv iHonda et al.l ( 2009f ) ob- 
servationally. The 3 /im scattered light comes from the 
layer to which the optical depth for the 3 /im radiation 
from the central star is about unity. Therefore, we define 
the layer as the disk 'surface' or '3 /im surface' in this 
paper. Next, we consider the ice-condensation front in 
the (R, Z)-plane. Without the accretion viscous heating, 
the temperature decreases along the radial coordinate 
R and increases along the vertical coordinate Z. Thus, 
the ice-condensation front becomes lower Z for smaller 
R and higher Z for larger R (see Fig. [3]). Then, there 
may be a cross-point between the 3 /im surface and the 
ice-condensation front. If we observe a disk with the 
scattered light, we will observe the 3 /im feature only 
from the outside of the cross-point. Therefore, we define 
the cross-point as the surface snow line in this paper. 
Namely, the surface snow line in this paper is the snow 
line which would be observed by the 3 /im scattered light. 

3. RESULT 

3.1. Structure of the ice-condensation front and the 
surface snow line 

Figure [3] shows the ice-condensation front (solid curve) 
and the surface snow line position (vertical tick) in a disk 
around a star with the effective temperature T* = 6000 
K, the luminosity L* = 10 Lq, and the mass = 1 
Mq. The layers of = 0.1, 1, and 10 are also shown 
as the dotted curves, where is the optical depth at 
the wavelength 3 /im from the central star. The layer 
of T 3Aim = 1 is the 3 /im surface of the disk. We also 
show the pressure scale-height calculated by the mid- 
plane temperature as the dashed curve in each panel. 

Comparing the panel (a), where photodesorption is 
omitted, with the panel (b), where the process works, we 
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Figure 2. Vertical profiles of disk annuli with the radius of (a) 10 AU and (b) 100 AU. The stellar effective temperature of 5000 K, 
luminosity of 10 Lq, and mass of 1 Mq are assumed. The top panels show the temperature profiles. The middle panels show the profiles 
of the changing rates of the grain size s due to condensation (double-dot-dashed), sublimation (dotted), and photodesorption (solid). The 
dashed and dot-dashed curves represent the rates due to photodesorption induced by the stellar (photospheric and excess) FUV luminosity 
and the interstellar FUV radiation, respectively. The bottom panels show the profile of the ice condensation fraction x. The horizontal 
axis represents altitudes from the disk mid-plane scaled by the radius of the annuli. 




R [AU] R [AU] 

Figure 3. The structure of the ice-condensation front in a circumstellar disk. The horizontal and vertical axes represent the distance from 
the central star R and the height from the mid-plane Z in astronomical unit, respectively. The panels (a) and (b) are the cases without 
and with photodesorption, respectively. The stellar parameters assumed are the effective temperature of 6000 K, the luminosity of 10 Lq 
and the mass of 1 Mq. The dust particle size is assumed to be 1.0 /un. The solid curve represents the condensation front below which 
water ice particles exist. The dashed curve represents the pressure scale height of the disk medium. The upper (thin), middle (thick), 
and lower (thin) dotted curves represent layers at which the optical depth at wavelength 3 lira from the central star becomes 0.1, 1 (3 /im 
'surface'), and 10, respectively. The vertical tick indicates the cross-point between the ice-condensation front and the 3-/im disk surface, 
i.e. the surface snow line position. 
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find that photodesorption wipes out water ice particles 
above the disk surface and the ice-condensation front al- 
most coincides with the surface throughout the disk. The 
position of the surface snow line (i.e. the cross-point of 
the ice-condensation front and the 3 /xm surface) shifts 
outward by a factor of about 4 for this case. As shown 
later (Figs. H]and[5|), the radial shift becomes larger and 
larger as the effective temperature increases. Therefore, 
photodesorption is very important to determine the po- 
sition of the surface snow line for Herbig stars. 

3.2. Stellar parameter dependence of the surface snow 

line 

In this section, we examine the dependence of the sur- 
face snow line position R sn ow on the stellar effective tem- 
perature T», the luminosity L„, and the mass M». In 
this calculation, we omit opacity of water ice particles for 
the computational cost. This simplification gives 20-30% 
smaller radius of the surface snow line position as shown 
in Figure |U This relatively small difference is probably 
caused by an offset of two opposite effects; (1) A smaller 
opacity in the case without ice opacity results in a higher 
temperature in the disk interior (i.e. a higher pressure 
height found in Fig. @| and also a higher surface height 
defined by T3 M , n = 1 (as also found in Fig. 0}. However, 
this leads a smaller density at the surface, and then, the 
condensation rate becomes slower and the surface snow 
line shifts outward. (2) The smaller opacity also results 
in a larger column density from infinity to the surface in 
order to account for r 3Atm = 1 at the surface. It leads a 
lower surface height and a larger density at the surface. 
Then, the condensation rate becomes faster, and the sur- 
face snow line shifts inward. This second effect exceeds 
the first effect in the cases shown in Figure |U In any 
case, the omission of the ice opacity does not change the 
mechanism determining the surface snow line position. 

3.2.1. Effective temperature 

As found in Figure HI the position of the surface snow 
line strongly depends on the stellar effective tempera- 
ture; as the temperature increases, the surface snow line 
shifts outward greatly and eventually disappears. Figure 
[5] shows this dependence more in detail for the cases with 
M, = 1 M and £* = 1 L© and 10 L Q . The dot-dashed 
curves indicate the cases without photodesorption, i.e., 
-Rsnow is determined by the balance between sublima- 
tion and condensation. On the other hand, the dotted 
curves are the cases without sublimation, i.e., R sn0 w is 
determined by the balance between photodesorption and 
condensation. The solid curves are the cases with all the 
three processes. The dashed curves are same as the solid 
curves but without the stellar excess FUV (see §2.4). 

For the sublimation-condensation cases (dot-dashed 
curves), i? S now is about 7 or 25 AU and almost inde- 

1/2 

pendent of T*, but depends on as R sn0 w oc LJ . 
This is because the gas temperature for the sublimation- 
condensation balance is almost constant at ~ 100 K at 
the disk surface where the dust temperature is deter- 
mined only by bolometric flux (i.e. oc L*/i? 2 ) and not 
the spectrum (see also §4.2). For the photodesorption- 
condensation cases (dotted curves), i? sn0 w strongly de- 
pends on as R sn0 w oc T} 2 . This strong dependence 
comes from the T* dependence in the FUV photon num- 



ber flux (see §4.2). It also depends on L» as -R sn0 w oc L 2 
which is different from the sublimation-condensation 
case. 

When all the three processes are taken into ac- 
count (solid curves), i? S now traces the sublimation- 
condensation case for lower T* and the photodesorption- 
condensation case for higher T*. In other 

words, the mechanism determining R sn0 w is the 
sublimation-condensation balance for lower T* 
and the photodesorption-condensation balance for 
higher T*. We here define a critical temperature 
T c as T* at the intersection point between the two 
curves of the sublimation-condensation case and the 
photodesorption-condensation case. If T* > T c , pho- 
todesorption is efficient and the surface snow line moves 
to the outermost part of the disk or disappears. Thus, 
T c is an indicator for judging the possibility to detect the 
3 /im water ice absorption feature in the disk scattered 
radiation (see §4.1). 

If we compare the dashed curves (the cases without 
the excess FUV) with the solid curves, we can see the 
effect of the stellar FUV excess on the surface snow line. 
Although it is still uncertain how strong the FUV excess 
is, we assume that the FUV excess accounts for only 
0.1% of L* (§2.4). We find an effect of the FUV excess 
in the L* = 10 L© case. If the FUV excess is stronger 
than this level, we should take it into account in order to 
determine the position of the surface snow line correctly. 

3.2.2. Stellar luminosity 

Figure [6] shows R sn ow as a function of T* with various 
Z/»s and A/* = 1 M©. It is seen that i? S now shifts outward 
as L* increases. This is because both the FUV radiation 
flux from the central star and the temperature in the 
disk increase as L, increases. It is also seen that the 
critical temperature T c becomes lower as L* increases. 

The dependence is approximately T c oc L* (see §4.2). 

3.2.3. Stellar mass 

Figure [7] shows R sn0 w as a function of T* with L* = 10 
Lq and various M*s. It is seen that the critical tempera- 
ture T c becomes higher as M* increases. This is because 
the condensation rate is approximately proportional to 
the square root of M„; the condensation rate depends 
linearly on the gas pressure with the assumption of uni- 
formly constant number fraction of water molecules in 
the whole disk, and the gas pressure depends approxi- 
mately linearly on the square root of M* (see eqs. A1-A3, 
and also §4.2). 

4. DISCUSSION 

4.1. Implications for the observation 

In this section, we discuss the effect of photodesorp- 
tio n on observ i ng the surface snow line. As proposed 
bv llnoue et all (|200l . we can resolve the surface snow 
line by using the 3 /Ltm water ice absorption feature in 
the scattered light if we have a good enough spatial res- 
olution. This method implicitly assumes the presence of 
water ice particles at the 3 /im disk surface. However, an 
efficient photodesorption may remove the particles from 
the surface and the absorption f eature from the scatt ered 
radiation. On the other hand, iHonda et al.l (|2009l ) de- 
tected the absorption feature in the scattered light from 
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Figure 4. Same as Figure \3\ but for various stellar effective temperatures. The thick and thin curves are the cases with and without the 
water ice opacity, respectively. Note that there is a difference in the axis scales between the panels (a) and (b) and the panels (c) and (d). 
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Figure 5. Position of the surface snow line as a function of stellar 
effective temperature with the stellar mass of 1 Mq . The horizon- 
tal and vertical axes represent the effective temperature and the 
position of the surface snow line, respectively. The thick and thin 
curves show different stellar luminosities as indicated in the panel. 
The dot-dashed curves show the snow line position where the equi- 
librium between sublimation and condensation is established. The 
dotted curves show the position where the equilibrium between 
photodesorption and condensation is established. The solid curves 
show the position determined by the equilibrium among photodes- 
orption, condensation, and sublimation. The dashed curves are the 
same as the solid curves but without the stellar UV excess. 



the outer disk at R > 140 AU of HD142527, although 
they did not resolve the snow line position because of in- 
sufficient spatial resolution. Therefore, photodesorption 
does not remove water ice particles completely from the 
disk surface of HD142527. 

As found in §3.2.1, photodesorption surpasses subli- 
mation at the disk surface and i? sn0 w is determined by 
the balance between photodesorption and condensation 
when T» exceeds a certain value which is defined as 
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Figure 6. Same as Figure [5] but for various stellar luminosities. 
The open circles indicate the critical temperatures T c defined in 
§3.2. 

the critical temperature T c . This temperature also de- 
pends on as found in §3.2.2. Let us over-plot the 
boundaries at which the mechanism determining R sn ow 
changes from the sublimation-condensation balance to 
the photodesorption-condensation balance on the HR di- 
agram and compare them with Herbig Ae/Be stars on 
the diagram (Fig. [3]). In particular, HD142527, where 
the water ice feature in the scattered light was detected, 
is the key object in this comparison. In addition, we also 
examine effects of other parameters such as photodesorp- 
tion yield, H2O abundance, etc. on the boundaries. 

As a reference, we define the fiducial model where 
Ah = 3 M©, the surface density at 1 AU Siau = 
1700 g cm -2 , the total (vapor+ice) number fraction of 
water X^® x = 1.2 x 10~ 3 , and the photodesorption yield 
Ypd = 1 x 10~ 3 . In addition, we consider three other 
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Figure 7. Same as Figure [5] but for various stellar masses. The 
stellar luminosity of 10 Lq is assumed. The solid, dashed, dotted, 
and dot-dashed curves show the surface snow line position with the 
stellar mass of 0.5 Mq, 1 Mq, 2 Mq, and 4 Mq, respectively. The 
open circles indicate the critical temperatures T c defined in §3.2. 

models: A, B, and C, where A t ^ 2 t ° is increased by a fac- 
tor of 2, Ypd is reduced by a factor of 2, and X^^ is 
increased by a factor of 2 and Y p d is reduced by a factor 
of 2, respectively (Table 1). In Figure we show the 
boundaries of these models by the dotted lines with the 
labels like 'fiducial' on the HR diagram. Interestingly, 
the boundaries of the models A and B are completely 
overlapped. The reason is the offset of twice X^ 2 t ® in 
the model A and half Ypd in the model B as shown in 
§4.2 analytically. In Figure [5J we also show the evolu- 
tion tracks of pr e-main sequence stars with various stel- 
lar masses from lYi et a l. (2001) (solid lines with labels 
of M*) and the stella r parameters of the obse rved Her- 
big Ae/Be stars from Ivan Boekel et ail (|2005f ) (circles). 
HD142527 is indicated by the square. 

In the upper-left region from the boundaries in Fig- 
ure [SJ the photodesorption-condensation balance deter- 
mines i? sri ow On the other hand, in the lower-right re- 
gion from the boundaries, the sublimation-condensation 
balance determines it. In the photodesorption- 
condensation case, the surface snow line moves out- 
ward significantly as shown in §3.2. Since iHonda et aLl 
(2009) detected ice particles on the disk surface at 
R > 140 AU of HD142527, we expect that HD142527 
is in the sublimation-condensation area. However, it is 
in the photodesorption-condensation area for the fidu- 
cial model. On the other hand, this object is in the 
sublimation-condensation area for the model C and 
marginally for the models A and B. 

Even if the surface snow line is controlled by photodes- 
orption, it is possible that R sn0 w < 140 AU as shown in 
Figures [5] and [6l Then, we show the boundaries at which 
R snow = 100 AU on the HR diagram in Figure El The 
solid lines with labels fiducial, A/B, and C are the bound- 
aries for the four models in Table [TJ We also show the 
boundary of 1000 AU for the fiducial model as the dashed 
line for a comparison. We find HD142527 between the 
boundaries of 100 AU for the fiducial and A/B models. 
In fact, this object is on the boundary of 300 AU for the 
fiducial model which is not shown in Figure |H1 Given 
the uncertainty of our simple model, we conclude that 
R snow in the HD142527 disk is 100- 300 AU which is sti ll 
consistent with the observations by IHonda et al.l (|2009t) . 
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Figure 8. Boundaries (dotted lines with fiducial, A/B, and C) 
at which the mechanism determining the position of the surface 
snow line changes. In the upper left side from the boundaries, 
the mechanism is the balance between photodesorption and con- 
densation. In the lower right side, the mechanism is the balance 
between sublimation and condensation. The labels fiducial, A/B, 
and C correspond to the models summarized in Table [JJ Note that 
the models A and B are completely overlapped. The inset shows 
the slope of the boundaries analytically derived in §4.2. The solid 
curves with labels of the stellar mass show the evolution tracks 
of pre-main sequence stars from Yi et al. (2001). The circles are 
the stellar parameters of the Herbig Ae/Be stars from van Boekel 
et al. (2005). The square shows HD 142527 observed by Honda et 
al. (2009). 

On the other hand, there are many Herbig Ae/Be stars 
in the upper-left region far from the boundaries of 100 
AU and even 1000 AU in Figure [9l Photodesorption is 
probably very efficient for these objects, and water ice 
particles disappear from the surface of the disk around 
the objects. Then, we expect that there is not the 3 /jm 
water ice feature in the scattered light from the disks. 
This will be an observational test to examine whether 
photodesorption is active or not on the disk surface. 

Finally we note a possible caveat in resolving the sur- 
face snow line by using the water ice 3 /im feature 
in the scattered light. As found in Figure HI the ice- 
condensation front and the 3 /Ltm surface (i.e. T3 Mm = 1 
layer) are always close to each other when T* < T c . When 
we are moving from the outer disk to the inner one on 
the scattered light image, thus, the strength of the 3 /xm 
feature may not disappear suddenly at the surface snow 
line, but may weaken gradually. In other words, the sur- 
face snow line is not a 'clear' line. In order to predict 
the change of the feature strength, we need to simulate 
observations, which would be a future work. 

4.2. Analytic expression of the critical effective 
temperature 

Here we try to derive an analytic formula describing 
the critical temperature, T c . The mechanism determin- 
ing i? sn0 w changes depending on T*; if T* < T c , it is the 
balance between condensation and sublimation, whereas 
if T* > T c , it is the balance between condensation and 
photodesorption. Thus, equation ([2]) can be approxi- 
mated to 

J Scon "i" Ssub ~ (T# < T^) 

\ Scon + Spd ~ (T* > T c ) ^ I 

When T* < T c , equation (0) becomes Xf^P w P sat , 
where X f ?* is the number fraction of H 2 molecules 
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Table 1 

Parameters for the (T C ,L») curves in Figures l8l and l9l 



Model 


Siau (g cm" 2 ) 


M, (M ) 


Y H 2 








fiducial model 


1700 


3 


1.2 x 10" 


a 


1.0 x 10 


-3 


A 


1700 


3 


2.4 x 10" 


-3 


1.0 x 10 


-3 


B 


1700 


3 


1.2 x 10" 


-3 


0.5 x 10 


-3 


C 


1700 


3 


2.4 x 10" 


-3 


0.5 x 10 


-3 
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Figure 9. Same as Figure[8] but the solid lines with fiducial, A/B, 
and C are the sets of the stellar luminosity and effective tempera- 
ture which give the position of the surface snow line -Rsnow = 100 
AU. The dashed line is the sets which give Rsnow = 1000 AU for 
the fiducial model. 

relative to all gas particles, P and P sat are the gas pres- 
sure and the saturation pressure of water vapor, respec- 
tively (see eqs. [B1-B5]). We have assumed that the 
vapor fraction X^° r ~ X^ aV The saturation pres- 
sure is a rapid function of gas temperature T. Thus, 
x f£\ p ~ p sat is established only when T ~ 100 K on 
the disk surface. This means that i? S now is the position 
where the dust temperature becomes T ~ 100 K. The 
dust temperature on the disk surface is determined by 
oT A w L^/AnR 2 . For a fixed T, we obtain Rs now deter- 
mined by the condensation-sublimation balance as 

-^snow,sub OC t J * ■ (9) 

This relation is very consistent with the dot-dashed lines 
in Figure [5] 

When T» > T c , equation JSJ becomes Ypd-Wabs oc 

IT Q 

^totai/ ' wriere V P d is the photodesorption yield, 7V a b s 
is the number flux of FUV photons absorbed by water 
ice grains, and p is the gas mass density (see eqs. [B4- 
B6]). The FUV number flux V a b s given by equation (B7) 
can be approximated to V a b s oc L^T^R for a fixed 
L» when the FUV absorption efficiency parameter (ra- 
tio of the absorption and geometrical cross sections of 
grains) Q a b s ~ 1 and the stellar FUV excess is omit- 
ted. The gas density p should be that at the 3 /jm sur- 
face: /93 Mm = p(T^^ m — 1). This can be approximated 

to p3[im °c oc M^ 2 R~ 3 ^ 2 , where Hp is the pres- 

sure scale-height. We have omitted the R dependence 
in the equatorial temperatu re Th because To is almos t 
constant at R ~ 100 AU (|Chiang fc Goldreichl fl997h . 
Note that /f^m is independent of the normalization of 
the surface mass density £ because we always see the 
layer with = 1. However, it depends on Hp which 
determines the density profile. As a result, we obtain 



Y pd L*T?R- 2 oc X^°Ml /2 R~ 3 / 2 , and then, R snow de- 
termined by the condensation-photodesorption balance 
becomes 

i?snow, P d cc T^LlY^Xf^M" 1 , (10) 
which agrees with the dotted lines in Figure [5] 

Finally, W6 define T c £lS when -Rsnow, sub — -Rsnow, pd- 

Then we obtain 

T G *L?'*Ml'^'*X*£'\ (11) 

Comparing this with the results shown in Figure [8j we 
obtain 

/ U V 1/8 / M. \ 1/12 
T c = 8500 K — — — 

"(t^) (t^J ■ (12) 

This formula would be useful for the readers to expect 
which balance determines R sn0 w If T* > T Cl one may 
expect the balance between condensation and photodes- 
orption. For the opposite case, the balance between con- 
densation and sublimation is expected. 

4.3. Effects of processes omitted 
4.3.1. Photodissociation 

We have omitted photodissociation of H 2 molecules, 
whereas this process is sugg ested to be important in 
disks around He rbig stars (jPontoppidan et ahl 120101 : 
IFedele et al.H201ll ). Here we justify this omission. 

The photodissociation rate is given by fc p dis ~ 
CpdisVpuv! where <7 p dis is the cross section for photodis- 
sociation and Np\]y is the number flux of FUV photons0 
According to Eel (|1984f ). a pdis ~ 10~ 18 cm 2 for H 2 
molecules. On the other hand, the photodesorption rate 
is given by fc pdc s ~ ^s 2 Q a , bs Y pd N FVV , where s is the 
grain size, Q a b s ~ 1 is the absorption efficiency param- 
eter, and Ypd is the photodesorption yield. Since the 
time-scale is r — 1/fc, the ratio of the two time-scales is 

Tpdis/Tpdcs = fcpdcs/fcpdis ~ 7rS 2 Ypd/cT p dis ~ 10 7 for s ~ 1 

pm and V p d ~ 10~ 3 . Then we find that photodissocia- 
tion is several orders of magnitude slower than photodes- 
orption. Therefore, photodissociation does not affect the 
balance between photodesorption and condensation. 

We may still argue the importance of photodissocia- 
tion in observations. When photodesorption becomes 
faster than condensation, water ice particles are fi- 
nally evaporated. This destruction time-scale is given 
by Tdcst = s/|s p dj, where s p d is the size changing 
rate by photodesorption (eq. B6). The ratio of this 

2 ^abs m the previous subsection is equal to Qabs^FUV- 
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time-scale to the photodissociation one is Tdcst/Tpdis ~ 
o"pdis/Oicos/(m H2 o^pd) ~ 10 5 for s ~ 1 /im and Y pd ~ 
10~ 3 (pice is the ice material density and toh 2 is the 
mass of a H 2 molecule). Thus, when the ice destruction 
is completed, evaporated H2O molecules have been al- 
ready photodissociated. Therefore, we expect that there 
is no H2O molecules in gas phase above the disk surface 
where photodesorption destroys ice particles. 

4.3.2. X-ray photodesorption 

We have negl ected desorption by X-ray (e.g., 
iNajita et al.l 120011 ) . The photodesorption rates by X- 
ray and FUV are given by fci « irs 2 QiYiNi, where i 
is X-ray or FUV, s is the grain size, Qi is the absorp- 
tion efficiency, Yi is the photodesorption yield, and Ni is 
the photon number flux. Thus, the ratio is fcx/fcFUV ~ 
(Qx/Qfuv)(Yx/Y fijv )(N x /N fvv ). 

X-ray luminosity of Herbig stars is obser ved in a range 
of 10 ~ 7 -10~ 4 of the bolometric luminosity (Stelz er et al.l 
2005). Given about two orders of magnitude higher en- 
ergy per a photon in X-ray than in FUV and that ~ 10% 
of the bolometric energy is emitted in FUV for Herbig 
stars, we obtain Nx/Npuy < 10~ 5 . The X-ray absorp- 
tion efficiency is Qx ~ 10 _1 which is estimated from the 
sum of atomic photoionization cross sections. Then, we 
find Qx/Qfvv ~ 10 _1 . The yield for X-r ay Yx is un- 
certai n for water ice, while Vpuv ~ 10~ 3 . INajita et al.l 
(|2001[ ) estimated Y x ~ 10-10 3 for CO from a consider- 
ation of 'spot heating' of subunits which is hit by a X- 
ray photon and is poorly connected thermally with other 
parts in a grain. They considered < 0.01 (im subgrains, 
which cover the parent grain, as the subunits. Compar- 
ing CO to H 2 0, the la tter has 5 times larger bindin g 
energy than the former (jSandford &: Allamand ola 1988), 
which reduces Yx for H2O by the same factor. In ad- 
dition, the connecting 'neck' between H2O ice particles 
is thic ker than that of silica considered in iNaiita et al.l 
(200l() because of sintering (|Siron 3 HOllt), further re- 
ducing Yx because of an enhanced conductivity of heat. 
Thus, we estimate Vx/Yfuv < 10 5 , although a more 
detailed study is encouraged. In summary, we obtain 
fcx/fcFUV < 10 _1 and we can omit X-ray photodesorp- 
tion at the disk surface. However, we note that X-ray 
photodesorption should be important in the disk inte- 
rior, where kx /fcFUV > 1 because iVpuv is significantly 
smaller than that at the surface. 

4.3.3. Dust settling and turbulent mixing 

In this paper, we have not considered dust settling and 
turbulent mixing although they may affect the vertical 
distribution of water ice particles and the surface snow 
line. If the dust settling toward the mid-plane occurs, 
the height of the 3 /im surface is reduced. This results in 
a larger density and a faster condensation at the surface. 
Thus, the surface snow line moves inward. To estimate 
the amount of the shift, we need further quantitative 
studies. 

The turbulent mixing may supply water ice particles 
from the disk interior to the disk surface even if photodes- 
orption removes the particles from the surface. However, 
the mixing time-scale is 1 x 10 3 yr at 10 AU and 1 x 10 4 
yr at 100 AU according to iHeinzeller et aTl ()201lD . while 
the destruction time-scale of 1 /jm ice particles by pho- 
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Figure 10. Position of the surface snow line as a function of dust 
grain size for the case of the effective temperature of 7000 K, the 
luminosity of 10 Lq, and the mass of 1 Mq. 
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Figure 11. (a) Height of the 3 /an surface relative to the 1 /un- 
sized grain case, (b) Height of the ice-condensation front relative 
to the 1 /xm-sized grain case. The dotted and dashed lines are the 
cases of 0.1 /tm-sized and 10 itm-sized grain cases, respectively. 

todesorption is 3 x 10 2 yr at 10 AU and 3 x 10 3 yr at 100 
AU. Therefore, the ice supply by turbulence is not fast 
enough. 

4.3.4. Grain size dependence 

Figure [Tul shows i? sn0 w as a function of dust grain size. 
We find that i? S now depends on the grain size in a com- 
plex way; the 1 /im-sized grain case gives the largest 
-Rsnow and smaller or larger grain size cases give smaller 
-RsnowS- This is because i? S now is the cross-point of the 
two curves, i.e. the ice-condensation front and the 3 /jm 
disk surface, which have a similar slope, and it is very 
sensitive to the relative position of the two curves. 

To see the reason of this complex behavior of i? S now > we 
show in Figure fTTI (a) the 3 /im surface height for the 0.1 
/im-sized grain and 10 /im-sized grain cases relative to 
the 1 /im-sized grain case. The panel (b) shows those of 
the ice-condensation front. From Figure ITT| we find that 
at around 10 AU, the 3 /im surface heights of both of the 
0.1 and 10 /im-sized cases are lower than that of the 1 
/im-sized case, while the ice-condensation heights of the 
both cases are similar to that of the 1 /im-sized case. This 
results in a smaller i? sn0 w for the 0.1 and 10 /im-sized case 
than for the 1 /im-sized case. This behavior is caused by a 
complex size dependence of the absorption and scattering 
opacities with and without water ice. More studies about 
this point are required in future. 

5. SUMMARY 

We have investigated the effect of photodesorption 
on the ice-condensation front and on the surface snow 
line position in a protoplanetary disk around a Her- 
big Ae/Be star. This study is inspired by the detec- 
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tion of an absorption feature of water ice at the wave- 
length 3 /im in the scatter ed light from the sur face of the 
disk around HD142527 bv lHonda et~aH (|2009l ). We have 
used 1+1D plane - parall el radiative trans fer alg orithm by 
IDullemond et al.l (|2002D and llnoue etail (|2009l ) to obtain 
the temperature and density structure in the disk. Then 
we have determined the amount of water ice particles by 
the balance between p hotodesorption, sublima tion, and 
condensation, following iGrigorieva et all f|2007f ) . 

We find that the effect of photodesorption depresses 
the ice-condensation front toward the mid-plane and 
pushes the surface snow line outward significantly. We 
also find that the surface snow line position is determined 
by the balance between photodesorption and condensa- 
tion for a higher effective temperature of the central star, 
while it is determined by the balance between sublima- 
tion and condensation for a lower effective temperature. 
We define the switching point of the two balances in the 
effective temperature as the critical temperature. The 
surface snow line is located at a small disk radius (< 
several tens of AU) if the stellar effective temperature 
is lower than the critical temperature. In the opposite 
case, the surface snow line moves outward rapidly and 
tends to be at a large disk radius (> several hundreds of 
AU) due to efficient photodesorption. We also derive an 
analytic formula to describe the critical temperature as 
a function of stellar and disk parameters. 



We discuss a few implications for the observation of the 
surface snow line in circumstellar disks. We plot the crit- 
ical temperature as a function of the stellar luminosity on 
the HR diagram. This indicates the boundary between 
two mechanisms determining the position of the surface 
snow line: the sublimation-condensation balance and the 
photodesorption-condensation balance. Most of Herbig 
Ae/Be stars are found in the photodesorption area. We 
also plot the line which gives the position of the surface 
snow line of 100 AU on the HR diagram. HD142527 is 
found slightly above the line, suggesting the surface snow 
line of the HD142527 disk is located at 100-300 AU from 
th e central star. This is consistent with the observations 
by iHonda et al.l (|2009T ) which indicated the surface snow 
line was located at < 140 AU. We also find that there 
are many Herbig Ae/Be stars whose snow line on the disk 
surface is > 100 AU or even > 1000 AU. It is difficult 
to detect water ice particles on the surface of the disk 
around these stars, which will be easily confirmed by ob- 
servations. Finally, we also find that the surface snow 
line position depends on the dust grain size by a com- 
plex way, which is worth investigating more thoroughly 
in future. 

We thank the referee for her/his constructive com- 
ments which were useful to improve the clarity of this 
manuscript very much. 



APPENDIX 
DISK STRUCTURE 

Here, we introduce details o f our calculations t o det e rmine the density a nd te mperature structure of a circu mstellar 
disk. Our method is based on IDullemond et al.1 (|2002fl . Ilnoue et al.1 (|2009f l. and lOka. Nakamoto fcldal (|2011h 



Hydrostatic equilibrium 

The dynamical structure of the disk in the vertical direction is assumed to be determined by the balance between the 
vertical component of the gravitational force by the central star and the gas pressure gradient (hydrostatic equilibrium) ; 

dP(R,Z) GM* 

—^— = - p{ R,Z)^ r Z, (Al) 

where P{R, Z), p{R, Z), G, and M* are the gas pressure, the gas volumetric density, the gravitational constant, and 
the mass of the central star, respectively. We adopt the equation of state of the ideal gas for the disk gas; 

P(R, Z) = Z ), (A2) 

where fc, T(i?, Z), p mc3 , n , and m u are the Boltzmann constant, the temperature of the disk gas, the mean molecular 
weight (we set this value to be 2.3), and the atomic mass unit, respectively. To obtain the vertical density profile 
p(R, Z) at a radius R by integrating equation (Al) with (A2) along the vertical axis Z ', we need the mid-plane density 
p(R, 0) as the normalization of p(R, Z). This is given by equating the surface density £(i?) with the integral of p(R, Z) 
along Z as 



S(-R) = / p{R,Z)dZ. (A3) 
The surface density S(i?) is given by equation (1). 



Radiative transfer 

The temperature in the disk is determined by the balance between the cooling and heating. The disk is heated by 
the irradiation from the central star. We ignore the heating by the viscous dissipation of the disk gas, because the 
snow line at the disk surface is located at relatively outer part of the disk where the viscous heating is negligibly small 
compared with the heating by the irradiation of the central star. The cooling rate in the disk is determined by the 
radiation energy transfer and the convection is not taken into account because of the low density in the disk. 



12 



Oka et al. 



Assuming a plane-parallel structure along the Z-axis in an annulus at a certain R, the equation of the radiative 
transfer is given by 

dZ 



A^J^ = P*u{S v {Z) - I,AZ)), (A4) 



where p, v, In,vi K>u> and S v are the cosine of the angle between the propagation direction and the Z-axis, the frequency, 
the specific intensity, the extinction coefficient, and the source function, respectively. The subscripts p and v of each 
quantity represent the direction cosine and the frequency of that quantity. If we assume an isotropic scattering and 
local thermodynamics equilibrium, the source function S v is given by 

S v {Z) = {l-m u {Z))B v [T)+m u {Z)^ J J^ + ^^ M (A5) 

where zv v , B V (T), and F^ 1 ^{Z) are the single scattering albedo, Planck function with temperature T, and radiation 
flux from the central star, respectively. 

We evaluate F- m U {Z) using the so-called grazing angle recipe (e.g., Chiang and Goldreich 1997, Dullemond et al. 
2002) as 

L v ( Tv (R, Z) \ 
W^ = 4^p«p^ —J, (Afl) 

where /3, L v , and t v (R, Z) are the grazing angle (i.e. incidence angle of the radiation from the central star into the 
disk surface), the luminosity of the central star, and the optical depth between a point (R,Z) and the infinity along 
the vertical direction (i?, oo), respectively. In this study, we set j3 = 0.05 throughout of the disk for simplicity which is 
a typical value of the grazing angle. We assume that the central star emits the black body radiation with an effective 
temperature T*. Hence, L v = An 2 R 2 B v (T*), where i?* is the radius of the photosphere of the central star. The optical 
depth t„ is defined as 

/•OO 

t v (R,Z) = / K uP (R,Z)dZ. (A7) 
J z 

The temperature of the dust particles are determined by the assumption of the radiative equilibrium; 

/•OO /> />oo 

girr + / pK v (l - vj v ) f I^ :U {Q.)dVLdv = 4?r / pn v {l - Wu)B u (T)dv, (A8) 
Jo J Jo 

where q m - is the heating rate by the irradiation of the central star and given by 

/•OO 

9irr = / p{Z)n v (l - m u )F- m ^(Z)dv. (A9) 



■Jo 

We can obtain the dif fuse radia tion field 1^.^, the temperature of the dust particles T, and the source function 5^(T) 
by solving equations (|A4[) and (|A8|) iteratively. Since a straight forward way to solve these equations converges very 
slowly, we employ a variable Eddington factor method (e.g., Dullemond et al. 2002; Inoue et al. 2009) for a rapid 
convergence. 

RATES FOR SUBLIMATION, CONDENSATION, AND PHOTODESORPTION 

Here, we describe the changing rate s of the radius of water i ce particles by sublimation, condensation, and photodcs- 
orption. These rates are according to IGrigorieva et al.l (|2007ft . 
The decreasing rate of the radius of water ice particles due to sublimation s su b is given by 

Ssuh = ~4^^ = (B1) 
$ sub = 4.08 x 10- 2 (^) (m 2 o^p\ g cm-V 1 , (B2) 



ltorr . 



where m S ub is the decreasing rate of the grain mass, s is the grain radius, p lcc is the density of water ice, $ su b is 
the mass sublimation rate, P sa t is the saturated water vapor pressure, r\ is the covering factor (fraction of the surface 
covered by sublimating material), and ^h 2 o is the atomic weight of the water molecule. We set p- lce to be 0.92 g cm" 3 
as Miyake & Nakagawa (1993 ) and consider that ther e exists purely icy dust particles (rj = 1). We adopt the saturated 
water vapor pressure same as IGrigorieva etaD (p00l : 



_ , 2.67 x 10 10 exp (- 6141 T 667K ) torr (T > 170K), 
sat ~ 1 5.69 x 10 12 exp (- 7043 - 51K ) torr (T < 170K). 



The increasing rate of the radius of water ice particles due to condensation is given by substituting the partial 
pressure of water vapor Ph 2 for the saturated water vapor pressure P sa t in equation (|B2I) . Namely, 



4npi ce s 2 Pice 



Photodesorption and Surface Snow Line in Protoplanetary Disks 
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$ con = 4.08 x 1(T 2 (m 2 o^j g cm^s- 1 , (B5) 

where m con and f coll are the increasing rate of the grain mass and the mass condensation rate onto the grains, 
respectively. The partial pressure of water vapor Ph 2 is given by the product of the number fraction of water vapor 
1h 2 o. vapor and the gas pressure P. 
The decreasing rate of the radius of water ice particles due to photodesorption is given by: 

•Vi = — (B6) 

^A'icc 

where toh 2 o = 3 x 10~ 23 g is the mass of a water molecule, Y" p d is the photodesorption yield, and -/V a b s is the number 
flux of the FUV radiation absorbed by ice particles. N a ^ s (R, Z) is given by 

N &hs (R, Z) = jf Qf° s dX, (B7) 

where F\(R, Z) is the incident FUV flux of the wavelength A at the position (R, Z) , an d Q^ bs is the absorption co e fficien t 
for water ice particles. We set the wavelength range of the integral of equation (IB7|) same as IGrigorieva et all {2007): 
X in = Q.Q91 /m i (13.6 eV) and A max = 0.24 /mi (5.1 eV). 

IWestlev et al.l (|1995al lbl) experimentally investigated the photodesorption of water ice and derived the yield as Y pd ~ 
1.0 x 10 -3 . lAndersson et al.l (|2006T ) and I Andersson fc van D ishocck (120081) carried out a cla ssical molecular dynamics 
simulation of photodesorption and derived the yield consistent with IWestlev et all (|1995al fbh . The weakness of their 
study was that their calculation was operated with a limited temperature and a limited wavelength range of the 
UV radiation. The latest experimental study on the photodesorption is lOberg etTa l. ( 2009) . They experimentally 
determined the rate of the photodesorption of water ice eliminating influences by other mechanisms and equipments, 
and obtained more detailed yield as a function of temperature and the thickness of the molecular layer of water ice as 
F p d(T, x > 8) = 10 _3 (1.3 + 0.032T) UV photon -1 . However, the results of these experimental and theoretical studies 
about photodesorption yield still have a large uncertainty. We set Ypd to be 1 x 10~ 3 as iDominik et al.1 (|2005[ ) and 
IGrigorieva et al.l (j2007h as the fiducial case. 
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